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Abstract 

We analyze the spontaneous U{V}r symmetry breaking at finite 
temperature for the simple O'Raifeartaigh-type model introduced in 
[1] in connection with spontaneous supersymmetry breaking. We cal- 
culate the finite temperature effective potential (free energy) to one 
loop order and study the thermal evolution of the model. We find that 
the R-symmetry breaking occurs through a second order phase tran- 
sition. Its associated meta-stable supersymmetry breaking vacuum is 
thermodynamically favored at high temperatures and the model re- 
mains trapped in this state by a potential barrier, as the temperature 
lowers all the way until T = 0. 



* Associated with CICBA 



1 Introduction 



It became clear after the work of Nelson and Seiberg [2J that global R- 
symmetry plays a key role in connection with supersymmetry breaking. In 
order to have spontaneous supersymmetry breaking at the ground state of 
generic models there must be a global U{1)r symmetry, but in order to have 
non-zero gaugino masses it is necessary that this symmetry be explicitly or 
spontaneously broken. The work of Intriligator, Seiberg, and Shih (ISS) |3] 
showed how this tension between i?-symmetry and supersymmetry can be 
exploited to find generic models with an acceptably long lived meta-stable 
supersymmetry breaking vacuum. Moreover, studying the Seiberg dual of 
M = 1 super-QCD it has been shown that, at high temperatures, the super- 
symmetry breaking vacua are dynamically favored over the "supersymmetry 
preserving" one^ so that the Universe would naturally have been driven into 
them [l]-[9], a possibility already discussed on general grounds a long time 
ago in [TU] . 

Different models with meta-stable symmetry breaking vacua and struc- 
tures rather different than those discussed by ISS have been also investigated, 
as for example those based in gauge mediation and extraordinary gauge me- 
diation, which cover a broad class of R-symmetric generic models with su- 
persymmetry breaking [TT]-|14j. 

There is a very practical mechanism proposed in [T] leading to sponta- 
neous U{1)r breaking. It applies to O'Raifeartaigh models with a continuous 
space of supersymmetry breaking vacua and degenerate tree-level vacuum 
energy. It has been shown in that work that, due to one loop corrections, 
spontaneous i?-symmetry breaking occurs a la Coleman- Weinberg in a very 
simple O'Raifeartaigh type model and for a wide range of parameters. More 
general models of this kind have been discussed in [15J and their thermal 
history has also been recently investigated ^16j . 

It is the purpose of this work to study the question of spontaneous U{1)r 
symmetry breaking at finite temperature and the resulting supersymmetry 
breaking pattern by analyzing the thermal evolution of the O'Raifeartaigh- 
type model introduced in [I]. To this end we will compute the finite tem- 
perature effective potential (i.e. the free energy density) by shifting as usual 
the relevant background fields and use the resulting quadratic terms (the 

^At finite temperature SUSY is always broken. With the quotation marks we mean the 
phase which, for zero temperature, corresponds to a supersymmetry preserving vacua. 
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mass terms) to perform the one-loop calculation. Studying numerically the 
corresponding one loop effective potential we will analyze the nature of the 
different phase transitions, showing how parameters of the model can be 
chosen so as to cover the desire range of critical temperatures at which R- 
symmetry breaking takes place. As we shall see, our numerical results are 
consistent with the general analysis presented in [TB] where a broad class 
of models for gauge mediation were considered. Indeed, in the classification 
of Extraordinary Gauge Mediation Models (EOGM) of [11], the model we 
analyze belongs to the type I class (provided one promotes the singlet mes- 
sengers to fields transforming in the 5 © 5 representation of SU{5)). Our 
analysis will confirm the thermal evolution scenario advanced in for type 

1 models, in particular concerning the existence of a metastable vacuum at 
high temperatures with no T = analog 

In the next section we introduce the model proposed in [Ij and describe 
its classical vacua, which includes a moduli space and a runaway direction. 
We then present the different terms that contribute to the one loop finite 
temperature effective potential V^jj. In section 3 we calculate V^jj along the 
pseudo-modulus, which is at the origin of the dynamically generated meta- 
stable vacuum, and analyze the R-symmetry breaking phase transition. We 
then extend in section 4 the calculation of V}jj by considering a background 
field that interpolates between the meta-stable vacuum and the runaway 
direction, and discuss in detail the resulting thermal scenario. We finally 
summarize and discuss our results in section 5. 

2 Set up of the model and the effective po- 
tential 

We consider the O'Raifeartaigh model for chiral superfields considered in pQ, 
with canonical Kalher potential and superpotential 

W = XX4>i4>2 + mi0i03 + ^m2(l)l + fX (1) 

This superpotential defines the underlying model which communicates super- 
symmetry breaking to the minimal supersymmetric Standard Model. Chiral 
superfields 0, {i = 1,2, 3) with R charges 

i?(0i) = -l, i?(02) = l, R{h) = s, (2) 
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represent the messengers of supersymmetry breaking and the spurion field 
X generates the model's pseudo-moduli space and has charge R{X) = 2. 
Parameters A, /, mi, and m2 will be taken, without loss of generality, as real 
positive numbers. 

The resulting scalar potential (we use the same notation for superfields 
and their lowest components) takes the form 

V''~'%X,<f),) = |A</.i02 + /P + |AX02 + mi03p+|AX0i + m2</)2p+|mi0i|2 (3) 
and its extrema consist of: 



a moduli space 
with 



dm) 



, arbitrary (4) 

V = f>0 (5) 



a runaway direction 

, (r) _ f \ ^ (r) _ ( fmiCps \ ^' Jr) 



~ \X'm,<pJ ' ^2 - V Am2 ^ ' 
A2/ 



oo, 



x« = ('4?#r (6) 



with 



V -^0. (7) 



The moduli space does not correspond to global minima of the potential but, 
as long as 

|X| < ^'-^ (8) 

where 

y = ^ (9) 

mim2 

it leads to local minima of the potential. Since the X field is i?-charged, such 
flat direction breaks the global U{1)r symmetry for any X ^ in the range 
(IE!). It is clear now that if quantum corrections produce a minimum at some 
point {X) 7^ of this flat direction, which then corresponds to a pseudo- 
moduli, the associated vacuum expectation value will spontaneously break 
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the R-symmetry. This was shown at T = in by computing the one-loop 
effective potentiaL We will now extend the analysis to include thermal effects 
by computing the finite temperature effective potential up to one loop, which 
takes the form [T7] 

V,'^iX'\ (Pf) = V'''%X''\ (t)f) + V^{X''\ (t)f) + V^{X'\ 0f ) . (10) 

The original fields are written in the form 

X = X'^ + x 

<P^ = <Pf + V^ (11) 

to proceed to compute the one-loop contribution by integrating terms quadratic 
in the fluctuations x, ipi. The zero temperature piece of the effective po- 
tential is given by the usual supersymmetric generalization of the Coleman- 
Weinberg formula 

^si^H-M^'og^ (12) 

where iSTr is the supertrace including a negative sign for fermions, Ai stands 
for the full mass matrix resulting from the shift ( fTTI) . M = M{X^'-'^,(f)f"-), 
and A is a mass scale. Concerning the finite temperature contribution, one 
has [nj 

= 1^ E f ds log (l T e-V^"^^) (13) 

where the sum is over all degrees of freedom {{rii} denotes the number of 
degrees of freedom, n = 2 for complex scalars and Weyl fermion and the upper 
(lower) sign is for bosons (fermions)). Finally, Aii denotes the eigenvalues of 
the A^-matrix. 

In order to make contact between the parameters of the model with scalar 
potential ([3]) and those of the Minimal Supersymmetric Standard Model 
(MSSM) one has to consider masses of the observable fields. It should be 
mentioned that a superpotential of the type ([H) should be in principle sup- 
plemented with a minimal gauge mediation (MGM) messenger 04 which, 
coupled to the spurion field X through a term of the form X(j)1, will ef- 
fectively give a mass to the otherwise massless gaugino [11]. Note that the 
introduction of this additional messenger would promote our model to a type 
III one, for which, instead of a condition of the form ([8]) stability requires 
an upper and a lower bound for X, X^ax > 1^1 > -^min, as noted in fTT] for 
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T = and discussed in [TB] for finite T. In the case of the model we consider 
one should adjust the parameters so that such bounds hold at all temper- 
atures and as the temperature grows Xniin(T) approaches the origin faster 
than the pseudomodulus minimum. We leave for a future work a detailed 
analysis of this issue and proceed to determine the orders of magnitude of 
the different superpotential parameters by analyzing sfermion masses. 

Sfermion masses can be extracted from the matter wave function 

renormalization through the formula |18j . 

where a is the running coupling constant of the underlying gauge theory 
(evaluated at the messenger scale, a/47r ~ 10"^) and 

i=l 

Here Aipi are the eigenvalues of the fermion mass matrix resulting from 
superpotential ([I]) and for simplicity we have set mi = m2 = m and defined 
X = XX /m. Given configuration @, Aip can be written in the form 

(XX* + 1 X 0\ 
X* XX* + 1 x\ (16) 
X* 1/ 

Formula f|T4l) is valid in the regime / -C for which supersymmetry is 
broken only in the effective field theory bellow the messenger scale by soft 
terms. 

Now one can check that 

N{x ^0) = X\ N{x^oo) = (17) 

Moreover, had we added the MGM messenger, the N behavior at infinity 
would have raised to so that we can take ~ in the whole range. 
In fact, if one scales X — > X/X and f Xf the coupling A completely 
disappears from superpotential ([1]) so that we can just set iV ~ 1 in f|T^ . 

Since one expects that the sfermion mass should be in the TeV scale, 
one infers from (HM that f/m 100 TeV, this in turn implying that m 3> 
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100 TeV. The estimate would remain nearly unchanged if instead of the 
assumption / ^ we consider the case / ~ m^. We conclude that for 
the analysis of the thermal evolution of the system, high temperatures will 
correspond to T ^ 100 TeV. 

3 The fate of the meta-stable vacuum 

We start by considering the effective potential for configuration (jl]), that is, 
we take (j)f = (f)f^^ = and X'^'' = X*^™) = X in formula (ITU]) . In this case 
the boson mass matrix takes the form (we omit the superscript m) 



Ml = 



while t 



j ml + X^X^ 


m2XX 








/A 


\ 


m2XX 


ml + X^X^ 


miXA 


/A 











mi AX 


2 

ml 














/A 





ml + X^X^ 


m2AX 





/A 








m2AX 


ml + X^X^ 


miAX 


\o 











miAX 


ml j 
(18) 


le fermion mass matrix reads 








j m{ + X^X^ 


m2AX 











\ 


m2XX 


ml + X^X^ 


mia;A 














mi AX 


m| 




















ml + A^X^ 


m2AX 














m2AX 


ml + A^X^ 


mi AX 













miAX 


ml j 










(19) 



Using this result, one can compute the zero-temperature one- loop contribu- 
tion (fT2!) . as originally calculated in [T], 

as well as the finite temperature one, eq. (fT3|) . which can be rewritten in the 
form 

V,^ = Ilj2 r ds (log(l - eV^^+^L/^^) _ iog(l + e-V^^+^k/^^)) 
^ i=i -^0 

(21) 
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One can scale X — rriiX/X and masses so that the effective potential 
only depends on the rescaled X and on two parameters: y, defined in eq. 
Qj, and r, given by 

r = ^ (22) 

mi 

so that Vi-^ = Vi^{X; r, y) with mi giving the mass scale. 

Eigenvalues M. Bi and M. pi (with i = 1, . . . , 6) of mass matrices ^A b and 
have to be computed numerically. Of course, at T = one reproduces 
the results in pQ thus finding that, for a wide range of parameters, there 
is a meta-stable vacuum where U{1)r is spontaneously broken. Concerning 
the thermal evolution we show in figure 1 the plot of V-^^ as a function 
of X for different temperatures. In figure 2 we represent the change with 
temperature of the region (shown in white) in the r, y plane where there is a 
U{1)r symmetry breaking local minimum of the potential satisfying (JHl). 




Figure 1: The effective potential as a function of |X| showing the second 
order phase transition (we have taken r = 4 and y = 0.2). The curve 
in the middle corresponds to the critical temperature which for the chosen 
parameters takes the value Tn/m = 0.95. 
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Figure 2: Plot oiy = \ f / {11111712) as a function of r = 1712/1711 for T = 0, 1, 1.5, 
and 1.8 (from left to right). The white region corresponds to a local (R- 
symmetry breaking) minimum (with no tachyons). 

Using different pairs of values (r, y) in the range where R symmetry break- 
ing occurs (white region in Figure 2) we have then found a second order phase 
transition at a certain critical temperature Tr, so that for T < Tr there is a 
minimum away from the origin, i.e. at X = {X) 7^ 0. 

Interestingly enough, changing parameters one can make the critical tem- 
perature vary in a wide range. For example, for the choice of parameters 
corresponding to Figure 1, {y = 0.2, r = 4) the critical temperature is 
Tr/iti = 0.95 while for y = 0.2, r = 2.07 it becomes Tfi/m ~ 10^^. In 
fact, by choosing parameters (r, y) closer and closer to the left frontier of the 
white region in Figure 2 one can lower the critical temperature as much as 
wanted. Taking into account the condition m 3> 100 TeV previously found 
from the requirement that m^/ ~ 1 TeV, we see that the critical temperature 
at which R-symmetry is broken can be adjusted in a wide range going for 
the two choices we have used as example, from Tr ^ 100 TeV to Tr ~ 1 
TeV. It should be noted that as the value of the critical temperature lowers 
the R-symmetry breaking VEV {X) gets closer to the origin. 

4 The fate of the runaway direction 

We will now study the behavior of the runaway direction as the temperature 
changes. To this end, we will follow an approach similar to that used in [1] 
in the case of the ISS and consider a path (X'"*, 0*"*) interpolating between 
the meta-stable supersymmetry vacua and the supersymmetric runaway di- 
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rection. A convenient choice of path is 



The function /i(03) should be chosen so as to concihate the behavior of X 
and 01 at the two-endpoints. An appropriate election is 

2 

h{y) = — a.Tcta.ncy (24) 

TT 

where c is a parameter to be chosen so that the path, which goes from the 
zero temperature meta-stable local minima (03 = 0) at X = {X) to the 
runaway value (03 — *• oo) does not have modes with negative square masses. 

We present in an Appendix the explicit form of boson and fermion masses 
for the path fl2^ . From their explicit form one can numerically study the 
effective potential as a function of 03 and the temperature, V^^ = V^^{(f)3, T), 
and determine the resulting minima landscape. First, one has to numerically 
compute the mass eigenvalues and then evaluate the zero temperature one- 
loop contribution to the effective potential (eq. (j20l) ) as well as the finite 
temperature one, V^, given by eq. (l2Ti) . 

One should note that at very high temperatures V^, as given by formula 
f fT3|) . becomes 



IT 
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V^' ~ -—T^ for T ^ oo (25) 
8 

Note that the negative sign in the effective potential is harmless since at finite 
temperature V^^ should be identified with the free energy as a function of 
the order parameter while the total energy is given by 

E = Kff-T^-^ (26) 

which is indeed positive for all temperatures. We show in figure 3 the free 
energy (left) and the total energy E (right) at very high temperatures. 
The figure clearly shows that although the energy is lower in what will become 
at zero temperature the runaway direction, the entropy contribution favors 
the non supersymmetric free energy minimum near the origin 
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From the numerical analysis of the complete effective potential ^^'^■^(^a, T) 
one infers the following scenario for the thermal evolution of the effective 
potential: 

• For T/m ^ 1 the potential has an absolute minimum at the origin in 
field space and it grows without bound for large values of 03. The zero- 

tcmpcraturc mcta-stable vacuum in the pscudomoduli direction has not 
yet started to develop and one finds, in addition, a local minimum at 
at a finite value 0^ (i.e. Vf*{(l)lT*) > Vf*{0,T*)) 

• As the temperature lowers, the slope of the potential at infinity de- 
creases until it becomes negative. The change of sign takes place at a 
temperature at which the absolute minimum of the potential is still 
at the origin. 

• At a lower temperature T5 the local minimum Vi"^* disappears. 

• At a lower temperature Tra, Vi^{(f)3 00, Tra) = V\^"^(0,Tra) so that 
the runaway minimum appears and a first order phase transition starts. 

• As already discussed, at a lower temperature Tr the i?-symmetry break- 
ing meta-stable vacuum arises. 

As an example, for the parameter choice r = 4,y = 0.2 already used to 
discuss the meta-stable vacuum evolution, the temperatures defined above 
take the values 

Th/m = 2.96 , Tb/m = 1.29 , Trjm = 1.14 , Tr/ui = 0.95 (27) 
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We have already described how changing parameters (r, y) towards the 
left border of the /^-symmetry braking region (white region in Fig. 2) lowers 
the critical temperature at which the transition to the meta-stable vacuum 
takes place. All other temperatures lower but their change is not so marked. 
As an example, for (r — 2.7, y — 0.2) one has 

Th/m = 1.5 , Tb/m = 0.99 , r^„/m = 0.81 , TR/m = 1 x 10"^ (28) 

Figure 4 shows a qualitative representation of the above scenario. 
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Figure 4: Evolution of the effective potential with temperature 
In order to exclude the possibility that the system escapes towards the 
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runaway direction instead of decaying into the meta-stable vacuum let us 
note that for T > Tf, the effective potential has an absolute minimum at the 
origin. Only for temperatures T < Ti, the runaway direction corresponds 
to an (asymptotic) global minimum of the effective potential. Since such 
temperatures are sufficiently low as to neglect thermal corrections, one can 
see PP that the barrier preventing the system to roll-down along the runaway 
direction has a width of order while its height is of order y^. Hence, 
by taking y sufficiently small the system will remain in the vacuum at the 
origin while Ti, < T < Tr and then smoothly evolve towards the meta-stable 
vacuum for T <Tr. 

5 Discussion 

We have analyzed the thermal evolution of the simplest O'Raifeartaigh-type 
model in which spontaneous R-symmetry breaking occurs dynamically, lead- 
ing to a runaway behavior at large fields and a meta-stable vacuum which, at 
zero temperature, spontaneously breaks supersymmetry. Studying the effec- 
tive potential at finite temperature we have shown that the f/(l)ij breaking 
arises through a second order phase transition. Remarkably, the critical tem- 
perature at which the R-symmetry breaking phase starts can be lowered by 
an appropriate choice of parameters and this also implies that the VEV of 
the spurion field X also decreases. 

We also analyzed the thermal evolution of the runaway direction find- 
ing, as expected, that high temperature contributions rise the asymptotic 
directions of the effective potential. Remarkably, we found that at high tem- 
peratures there is an extra local minimum of the effective potential, though 
energetically unfavored with respect to the meta-stable vacuum. At some 
temperature {Tb) this local minimum disappears. 

The whole thermal evolution sequence is as follows: At high tempera- 
tures the model is driven to the meta-stable SUSY-breaking vacuum. As the 
temperature decreases, the SUSY runaway direction becomes energetically 
favored but the transition between phases is long lived, so the system remains 
in the meta-stable vacuum. There is also an extra local minimum but with 
higher effective potential than the meta-stable vacuum. As the temperature 
decreases this extra minimum fades away. Finally, at an even lower temper- 
ature (Tr), the R symmetry is broken and a second-order phase transition 
occurs. This sequence, with the exception of the existence and eventual dis- 
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appearance of the extra local minimum, is similar to the one described in 
[1]-[S] for the magnetic dual of SuperQCD. As stated in the introduction, 
the model studied here can be extended to the form of a type I model in the 
classification of ref . [11] . The general properties of the thermal evolution of 
these models was discussed in [16] and our numerical analysis of the vacuum 
structure at different temperatures is consistent with them. In particular our 
results confirm the existence of an extra vacuum at high T in addition to the 
one at the origin, with no analog at T = 0. This extra vacuum disappears 
as the temperature lowers below Tb. 

An implicit assumption necessary to apply our results in a cosmological 
context is that the reheat temperature Treheat is larger enough (with respect to 
the supersymmetry breaking scale) as to guarantee that the supersymmetry 
breaking history develops quasi-statically, in a situation of thermal equilib- 
rium. This justifies to look for the minima of the free energy not taking into 
account possible interaction between fields and the heat bath. Ignoring the 
possibihty of non-equilibrium situations our results suggest that although the 
runaway direction starts to develop before the R-symmetry breaking meta- 
stable minimum appears, the system will not roll-down from the minimum 
at the origin because of the existence of a very high barrier so that when the 
R- symmetry breaking meta-stable vacuum is available, it will evolve to it 
and remain there for a sufficiently large time as to ensure that the Universe 
is still trapped there. 

We would like to end this work by pointing out two directions in which we 
hope to continue our investigation on R-symmetry breaking and supersym- 
metry breaking at finite temperature. One concerns the analysis of models 
with explicit R-symmetry breaking which, under certain conditions, have 
supersymmetric vacua, runaway directions and meta-stable vacua [IH]. As 
discussed in [20], the way in which R-symmetry is broken (spontaneously or 
explicitly) leaves a clear imprint on the phenomenology of the MSSM and it 
is then worthwhile to study broad classes of such models so as to compare the 
resulting thermal patterns. The other direction is related to the analysis in 
[2T| on how pseudomoduli arising in generalized O'Raifeartaigh models from 
additional global symmetries can be candidates to dark matter (see also [22]). 
In this context it would be of interest to investigate the thermal evolution of 
such models along the lines developed here. We hope to analyze these issues 
in a future work. 
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6 Appendix 

We write the boson and fermion mass matrices corresponding to the path 
(!23|l in the form 



■^^=1 A 



KAl I A B 

^-=i B A 



(29) 



where A and B are symmetric 4x4 matrices with nonzero elements 

_ r'y^/'h{<Ps?mf .,2/3,2/3 ,2/3„.4/3 

^11 - ^2/3^2/3 +y '^^ ^ ^1 



A 



12 



A 



21 



Ai4 = Ail = -m^'^^^v^v^ 



^2^2/3^^7/3 



22 



m? + y2/3^f A^/^mf + xoA - xo/i(03)A + 



/mT^fA^Y 



-423 = A32 = mir I xqA + -^^^ — {niiy + XQ\)h{(j)^) 



^34 


= ^43 


A44 


= ml 


B12 


= B21 


Bl3 


= -B31 


B23 


= -B32 



4/3 



+ XoXnii — xoA/i(03)mi 



/,2/3 



531 = mf xo^^A^/3(/i(03) - 1^ 
^32 = -mlry{h{(l)s) - 1) 
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